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Quantum mechanics in the general quantum systems (V): Hamiltonian eigenvalues 

Zhou Li and An Min Wan£| 

Quantum Theory Group, Department of Modem Physics, 
University of Science and Technology of China, Hefei, 230026, P.R.China^ 

We derive out a complete series expression of Hamiltonian eigenvalues without any approximation 
and cut in the general quantum systems based on Wang's formal framework [^. In particular, we 
then propose a calculating approach of eigenvalues of arbitrary Hamiltonian via solving an algebra 
equation satisfied by a kernal function, which involves the contributions from all order perturbations. 
In order to verify the validity of our expressions and reveal the power of our approach, we calculate 
the ground state energy of a quartic anharmonic oscillator and have obtained good enough results 
comparing with the known one. 

PACS numbers: 03.65.-w, 03.65.Ca 



I. INTRODUCTION 



To determine the Hamiltonian eigenvalues is a basic problem in quantum mechanics. As is well known, the exact 
solution of the Schrodinger equation can be obtained only in some special cases, i.e., for several elementary systems like 
• the hydrogen atom, the H2 molecule, or the harmonic oscillator. In the majority of cases, approximation techniques 
Oh' have to be employed in calculation of the Hamiltonian eigenvalues in the general quantum systems, and the precision 

> 1 and computability become focus. 
^ ^ Perturbation theory is one of the few principal methods of approximating solutions to eigenvalue problems in 
^ I quantum mechanics. The formalism of Rayleigh-Schrodinger perturbation expansion expresses an eigenvalue as a 
formal power series of the coupling constant A: 

00 

^total ^ ^ 

*^ ' m— 

■rj" At the 4th order perturbation level, the explicit form of C4 has appeared a little complicated. Nevertheless, for a 
' given nth order perturbation, one can obtain, in principle, the form of c„ according to Kato's [H or Bloch's 3 formal 
expression. 

The normal perturbation theory frequently meets two problems: one is about whether the power series is convergent 
in some neighborhood of A = or it is only asymptotic as A ^ 0; another is about the strong coupling regime. While 
On the problems of strong coupling is usually overcome with various kinds of renormalization techniques, for example, 
' the renormalization scheme recently used by Cizek and Vrscay 0, Q , summation techniques are employed to give a 
^ . divergent perturbation series any meaning beyond a mere formal expansion fsj. 

In our point of view, these problems closely connects with the whole and deep knowledge about the expansion ([T]) 
^ . , because the perturbation series is inherently multiple in the general quantum systems. In mathematics, a reasonable 
' rearrangement of a multiple series is often significant when a cut approximation needs to be introduced. Consequently, 
we wonder whether the expansion as a power series of the coupling constant is a unique choice, and whether a more 
explicit expression of total Hamiltonian eigenvalues exists. Moreover, we would like to find a systematical and new 
approach for the calculation of Hamiltonian eigenvalues in the general quantum systems. 

More than three year ago. An Min Wang, one of authors in this article, presented his research on a formal framework 
of quantum mechanics in the general quantum systems [H 0i 9 and made a conjecture about the total Hamiltonian 
eigenvalues. Just based on Wang's works, a complete series expression of Hamiltonian eigenvalues in a general quantum 
system without any approximation and cut is explicitly obtained by using some skills in mathematics and physics. 
This expression is simply not a power series of perturbed parameter, but a series of power of a kernal function as well as 
its derivatives that involves the contributions from all order perturbation. No cut and approximation are introduced, 
and the general term is given out. In special, our kernal function as well as its derivative is equal to O(A^). It implies 
that our series is obviously improved in its approximation content not only involving the higher order contributions 
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but also being more suitable to cut upto some given nth term since its decrease is more rapidly than the series in the 
normal theory. Of course, if expanding our expression of the total Hamiltonian according to the perturbed parameter, 
our expression is consistent with one in the normal perturbation theory. However, our expression really shows the 
physics nature because it directly comes from the law of quantum dynamics, but not a transcendent input ~ expansion 
according to the power of perturbed parameter which covers up a fact that perturbation series is inherently multiple 
in the general quantum systems. Moreover, it will be seen that the new conclusion can be obtained by our expression. 

It is worthy to point out that it is inevitable and acceptable that the complete expression of Hamiltonian eigenvalues, 
just like ours, is an infinite series, because that a general quantum system has usually no a compact solution without 
any approximation. More importantly and interestingly, from our expression we further present a practical approach 
to find the total Hamiltonian eigenvalues only through solving an algebra equation satisfied by our kernal function, 
which is specially suitable to calculation in a computer. In order to verify the validity of expression of Hamiltonian 
eigenvalues and check the accuracy of calculation of Hamiltonian eigenvalues using our approach, we study, as an 
example, a quartic anharmonic oscillator in the normalization, which is frequently used as a test-stone of new methods 
of calculating the eigenvalues. 

This article is organized in this way: in section |ll] we give an overview on how to get our expression of the total 
Hamiltonian eigenvalues in the general quantum systems. The detail is put in four appendixes; in secton lTTTl we present 
how to obtain the total Hamiltonian eigenvalues by solving a given algebraic equation; in section IIVI we calculate the 
ground state energy of the quartic anharmonic oscillator for different coupling constants; in section |Vl we make the 
summary. 



II. DERIVATION: HAMILTONIAN EIGENVALUES 



In the section, we would like to derive out an explicit expression of Hamiltonian eigenvalues in a general quantum 
system. Actually, our purpose is just to prove Wang's conjecture proposed in his theoretical framework of quantum 
mechanics in the general quantum systems [ij. 

In Wang's work, the complete series expression of time evolution operator or transition amplitude in a solvable 
representation {Hq]^'^) = Ej\^^)) is presented as 



where 
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It is clear that Wang's expression can be thought of to be exact in the sense that this series involves contributions from 
all order perturbations and has no any approximation. In above expression, for simplicity, we have not considered the 
degenerate cases. It is important that here we have used the following subtle method of dividing Hamiltonian matrix 
in the solvable Hq representation 



{H} = {Ho} + {H,} 
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Such an expression has separated the total Hamiltonian matrix into diagonal part and off-diagonal part, rather than 
nonperturbative part and perturbative part, which is marked with the coupling constant in the normal perturbation 
theory. It is never trivial in our derivations. Breaking the accustomed mentality in the normal perturbation theory 
and using the more formalized expansion form of time evolution operator are advantages of Wang's theory, and they 
bring us successfully to arrive at our purpose. 
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Obviously, there are many apparent singular points in the expression ^ of AJ'' , but they are fake in fact. In 
Wang's paper [l[ this problem has been fixed by finding their limitations in terms of contraction and anti-contraction 
of energy summation indexes. Here, we will theorize Wang's method and further prove Wang's conjecture about the 
eigenvalues of Hamiltonian. 

A key trick using here is that we start from the partition function and rewrite it by using Wang's expression ([2]) 



-iE-ft 



(6) 



where E-y are just the total Hamiltonian eigenvalues that wc would like to find. It is largely helpful for removing the 
unexpected "fake" singular points in Wang's framework, which can be seen in AppendixE) Another key skill is that 
we use the contraction and anti-contraction of energy summation indexes developed in [l| . It makes all the apparent 
singular points in the partition function expressed by Wang's framework are neatly removed, which can be seen in the 
Appendix[Bl It is interesting that, based on the proof in Appendix [Cl bv particularly analyzing, skillfully recombining 
the summations over energy indexes and perturbed order indexes, we arrive at 
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R-y{z) is a kernal function that involves contributions from all order perturbations, and plays a key role in our 
expression. It is clear that i?™+^(z) and d™i?™+^(z)/(iz"'|z=o have the same order of magnitude of perturbed 
parameter A, but and (d'"+ii?™+2(z)/dz"+i|^=o) / (d'"i?™+i(z)/dz"|^=o) is equal to 0{X^). This 

implies that the approximation ability of this series is obviously improved and then it is more suitable to cut upto 
some given the TOth term since it decreases more rapidly with A. 

However, more importantly and interestingly, we obtain 
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It is proved in the Appendix (|D|) by expending the partition function in Eq.(I7]) into the time power series and 
verifying the coefficient power relation. In fact, the form of Eq.® has its physics origin rather than the mathematics 
arbitrariness, and it is valid in the general quantum systems independent of the form of Hamiltonian. Therefore, we 
think that the complete series expression of Hamiltonian eigenvalues in the general quantum systems as below 
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Obviously, it is simply not a summation over the purturbed parameter, but a series of power of the kernal function 
R-y{z) as well as its derivative at z = 0. It is completely different from the normal perturbation theory in its 
thoughtway. In particular, when a cut is introduced in a practical calculation, a higher A^ term than the last term 
of cut part is dropped, but not a higher A term than the last term of cut part is dropped in the normal perturbation 
theory. 

It is easy to verify that our result is consistent with one in the normal perturbation theory if one expands our 
expression (|10p according to the order of perturbed parameter for the known lower order forms in the textbook and 
refs. However, our expression involves the contribution from all order perturbation and has a neat form of general 
term. Rearranging summation in our expression is helpful for theoretical derivation and practical calculation since its 
completeness, orderliness and clearness. In particular, it provides a physical reason how to chose part contributions 
from higher order perturbations, which is able to simplify the calculation and lead the result more precision. An 
example has been presented in Ref.[8]. In fact, our rearrangement summation is reasonable because it reveals the 
mathematical beauty and then physical nature. Moreover, its new content, at least, to surprise us, will be obtained 
in the following section. 

Of course, our expression of the total Hamiltonian eigenvalues can be thought of to be exact in the sense that this 
series involves the contributions from all order perturbations and has no any approximation and cut in form, as well 
as the general term is obtained. 
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III. CALCULATION: AN ALGEBRA EQUATION 



Although the complete series expression of total Hamiltonian eigenvalues (|10[) in the general quantum system 
arrives at our theoretical aspiration and we believe that it is interesting and important in the formulism of quantum 
mechanics, we have to admit that this expression is probably too complicated to be practical. It is inevitable and 
acceptable that the complete expression of Hamiltonian eigenvalues (|10p is an infinite series, because that a general 
quantum system has usually no a compact solution without any approximation. It seems to be not convenient in 
the calculation the nth order derivative of i?"^^, but because that Rj(z) is such a function with a product form 
— z)^ such a difficulty is not serious. However, except for a reasonable rearrangement summation and a neat 
general term to reasonably involve the higher order contributions, what is more in our expression than one in the 
normal perturbation theory for the calculation of total Hamiltonian eigenvalues. In this section, our purpose is just 
designed to answer this problem. 

Actually, we find that the difference between 7?'s and H^s eigenvalues IS.E^ is equal to the coefficient of 2;'s 0th 
power term in a Laurent series of F^[z) as following 



= zln 1 + 



R-y{z) 



(11) 



It is easy to be proved. In fact, setting u — R^{z)/ z, we can rewrite F^{z) = R-y{z) ln(l + u)/u. Obviously, the limit 
of ln(l + u) ju when u ^ is 1. Thus it means that ln(l + 7i)/u can be expanded as a Taylor's series at u = 0, that is 



F^{z) = R^{z) 



oo 



(-1)'= 



^ k 

k=0 



(12) 



Again note that i?^(z)'^+^ can be expanded as a Taylor's series at z = 0, we finish the proof of our above conclusion. 
Now our task is to seek an explicit form of Laurent series of Fj{z). We first set is a solution of following equation 



R-y{aj) 



(13) 



When a~f is finite, from above equation it follows that R-y{z) can be expended as 



Rj{z) = Rj{aj) + 
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Then, we rewrite 
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Since R-y(z)/{z — a-y) is canonical a,t z = a-y from Eq. (|13p for a finite a^, the second term in above equation does 
not involve z's 0th power part. Actually, this is a reason why to set Eq. (fT5| . While the first term in Eq.ljTS]) can be 
expanded as 



zln 



k=l 



Therefore the Laurent series of F^{z) has its coefficient of z's 0th power to be —a.y. It implies that we obtain 



A_E/-y E-y 



(16) 



(17) 



This conclusion is so interesting because it tells us that the total Hamiltonian eigenvalues in the general quantum 
system can be calculated through solving an algebra equation 



R^{-AE^) = AE^ 



(18) 



and then adding its solution AE^ to the Hq^s eigenvalues E^. More obviously, the total Hamiltonian eigenvalues E^ 
is a solution of following algebra equation 

E E - — ^..^ J..^ 1 . — ^ — ^ = 1 (19) 
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This conclusion distinctly reveals the advantages of Wang's formal framework and goes to an extent that we have 
ever not researched using the known other theory. It must be emphasized that the contributions from higher order 
even all order perturbation are naturally involved in such an algebra equation. The difference between our approach 
and usual one has been clearly seen here because that we give up the accustomed calculation method order by order 
in the normal perturbation theory. 



IV. AN EXAMPLE: QUARTIC ANHARMONIC OSCILLATOR 

In this section, we attempt to verify the validity of expression of Hamiltonian eigenvalues and check the accuracy 
of calculation of Hamiltonian eigenvalues using our method. As an example, we study a quartic anharmonic oscillator 
in the normalization, which is frequently used as a test-stone of new methods of calculating the eigenvalues. 

Harmonic oscillators and their anharmonic counter parts are extremely important model systems in all branches of 
quantum physics and particularly in quantum field theory. Even order anharmonic oscillators defined by 

=p2^a;2 + Ax2™,m = 2,3,4- •■ (20) 

were studied in Bender and Wu's seminal work [ol. [loL . and then Simon made a rigorous analysis of the mathematical 



property in |12l |. The perturbation expansions of the anharmonic oscillators diverge strongly for the coefficients c„ 
in Eq.([T]) always grow factorially. Therefore anharmonic oscillators are frequently used to test new approximation 
techniques. 

Let us begin to consider a quartic anharmonic oscillator. Its Hamiltonian reads 



da;^ 



^ = -332+^' + ^^' (21) 



Obviously, the primary task is to write out our kcrnal function R^. But, as it is defined, is the summation of 
an infinite series. Thus, a certain cut approximation should be made, that is, some terms of this series need to be 
omitted. So, what strategy should be adopted when trying to distinguish those of importance from trivial terms? 
One may think of the tactics when making perturbation approximation, we just pick those terms of low powers of the 
coupling constant, and drop the higher powers. We are going to point out that such a tactic won't be very proper 
here, for the two reasons below. First, the usual perturbation expansion is a power series of the coupling constant, 
but, our expression has no longer been simply a power series of coupling constant. Second, even though we admitted 
that the longer the term the smaller the value, the the longer the term, the bigger the number of such terms. So, 
we propose a new one, according to the fact that the states nearby exert larger impact than the states faraway in our 
expression, we drop the terms reflecting the effect of the states far away from the initial one. 

To put our idea into reality, we choose such some terms that the summations over 71 , 72 , • ' ■ j 7/ has a maximum 
value n and define them as a new series 

^ V V g77lg7l72...g7,7 

max{7i,72---7i}=n 

Thus, we can rewrite the kernal function Rj as the summation of R'jiz, n) over n as below 

00 

i?^(z) = ^i?^(z,n) (23) 

n=l 

In the calculation of the ground state energy of quartic anharmonc oscillator we make the cut approximation to N 

N 

i?^(z)«5]i?!;(z,n) (24) 

n=l 

To justify our approximation scheme, we've made some numerical calculation as shown in Table. [D And our result 
of the ground state energy of quartic anharmonic oscillor for different coupling constants is shown in Table. [H] 
These results are good enough comparing with the known one in [l^ . 
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TABLE L The value of 7?o(0,n). 



A 


n = 3 


n = 5 


n = 11 


n = 21 


n = 51 


n = 101 


n = 199 


n 1 
U.l 


-y.loo / e-UUo 


-4.0/ /oe-UU4 


-i.44yue-uu / 


-o.UiU / e-Uiz 


-o.iizye-uz4 


-z.oyiie-u4i 


-4.o4DOe-UDo 


0.2 


-3.10346-002 


-3.5160e-004 


-9.8794e-008 


-3.0709e-010 


-7.5163e-020 


-6.55726-034 


-4.48076-057 


1.0 


-3.4615e-001 


-2.1281e-002 


-8.6467e-004 


-1.5403e-006 


-6.5085e-013 


-4.51526-022 


-3.47686-040 


2.0 


-8.1818e-001 


-1.2014e-001 


-1.5617e-003 


-8.4901e-007 


-2.3659e-011 


-4.23706-018 


-4.47326-032 


10.0 


-4.7872e-H000 


-1.50666+000 


-2.0541e-001 


-5.5894e-004 


-2.8380e-006 


-4.80536-012 


-1.26836-020 


20.0 


-9.78266+000 


-3.42216+000 


-1.58296+000 


-7.2278e-001 


-1.3054e-003 


-3.13656-008 


-1.82636-014 


100.0 


-4.97796+001 


-1.89996+001 


-3.08936+001 


-1.29126+001 


-4.27346+002 


-3.67836-002 


-5.21846-006 



TABLE II: The value of the ground state energy of quartic anharmonic oscillator, using the the summation of the first n Ro{z, i) 
(i = 1, 2 ■ ■ ■ , n) as the approximation of R{){z). 



A 


n 


Eo 


A 


n 


Eo 


0.1 


10 


1.06528570130099 


0.2 


10 


1.11829330436519 




20 


1.06528550957781 




20 


1.11829265486895 




30 


1.06528550957275 




30 


1.11829265444366 




50 


1.06528550957275 




50 


1.11829265444348 




100 


1.06528550957275 




100 


1.11829265444348 




200 


1.06528550957275 




200 


1.11829265444348 




known" 


1.0652855095437176888 




known" 


1.1182926543670391534 


1.0 


10 


1.39337105560387 


2.0 


10 


1.61122760597946 




20 


1.39235392111137 




20 


1.60754799112121 




30 


1.39235164865408 




30 


1.60754155853087 




50 


1.39235164313030 




50 


1.60754130410410 




100 


1.39235164312960 




100 


1.60754130407997 




200 


1.39235164312960 




200 


1.60754130407997 




known" 


1.3923516415302918557 




known" 


1.6075413024685475387 


10.0 


10 


2.47630097947871 


20.0 


10 


3.18161125567721 




20 


2.45355539526673 




20 


3.02112722285399 




30 


2.44923642985496 




30 


3.01172336279951 




50 


2.44917490466071 




50 


3.00996284534114 




100 


2.44917407783815 




100 


3.00994481629290 




200 


2.44917407782312 




200 


3.00994481558327 




known" 


2.4491740721183869183 




known" 


3.0099448155577821983 


100.0 


10 

20 

30 

50 
100 
200 
known" 


8.08464496277487 
5.14809927717057 
5.02007347355405 
5.00376751877937 
4.99942534973870 
4.99941754801155 
4.9994175451375878293 









"Here, the value denoted by "known" from Ref. |13il . 



V. SUMMARY 

In this article, we proposed a new approach calculating the Hamiltonian eigenvalues in the general quantum systems. 
In theoretical form, we obtain a complete series expression of total Hamiltonian eigenvalues without any approximation 
and cut. In practical calculation, we present an algebra equation satisfied by the total Hamiltonian eigenvalues. These 
conclusions is based on Wang's framework of quantum mechanics in the general quantum system. Consequently we 
can say that Wang's formulism of quantum mechanics in the general quantum systems is useful. In fact, the revised 
Fermi's gold rule as well as its calculation Q also accounts for Wang's theory has itself advantages. 

Although only the non-degeneracy and discrete case is considered here, but our derivation can be extended to the 
continuous and/or degeneracy case. 

It must be emphasized that an accustomed mentality that expands the series according to the perturbed parameter 
in the normal perturbation theory is given up, the contributions from all order perturbation is involved via a kernal 
function, which plays a key role in our theoretical expression and our calculation approach. The advantages of 
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our kernal function have been mentioned in the previous sections. The calculation about this kernal function is an 
important task. However, it only involves a product of matrices, and then it can be efficiently calculated by using 
a computer. So in the computability and usability, our approach is not weaker. In formal beauty, please see our 
concrete expressions, no more words need to say. 

Perhaps, our expression is thought of a rearrangement of perturbation series. However, such a rearrangement 
summation in form is never trivial, it really shows the physics nature because it directly comes from the law of 
quantum dynamics, but not a transcendent input - expansion according to the power of perturbed parameter. In 
mathematics, the reasonable rearrangement of a multiply series is often significant if a cut approximation is needed. 
We think that the rearrangement summation in our expression make us reasonably involve the contributions from 
part higher even all order perturbations from our derivation. In fact, our expression is more neat, more explicit and 
more deep than the known Kato's and Bloch's Q. This should be a reason why it can lead to an interesting 
conclusion - the total Hamiltonian eigenvalues can be calculated by solving an algebra equation. 

Our approach to calculate the total Hamiltonian eigenvalues also gives up the accustomed way order by order, it 
makes the calculation to be simplified, but still involves the contributions from important part higher order pertur- 
bations when a cut approximation is introduced. Actually, it is similar to choose the important parts in higher order 
(equal to or more than the fourth order) contributions based on some physics reasons. Specially, our equation about 
the total Hamiltonian eigenvalues can more naturally and conveniently involve the contributions from higher even 
all order perturbations, and it is easy to be solved numerically in a computer. Therefore, our approach is able to 
remarkably simplify the calculation as well as advance the precision. 

As an example, we applied our approach in the calculation of the ground state energy of a quartic anharmonic 
oscillator. The highly accurate results for the energies of the ground state energy of quartic anharmonic oscillators of 
different coupling ranging from 0.1 to 100 are yielded. Something we should mention is that the factorial divergence 
in the usual perturbation expansion of even order anharmonic oscillators no longer appears in our expansion, and our 
way of "summing" expansion, which should be owned to the fact that our expression is complete, is rather simple and 
neat in form and also proved to be effective. On the basis of the results presented in this article, one can expect that 
our new method should also give good enough results in some other quantum mechanical problems, and we'll proceed 
our work in the near future. 
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APPENDIX A 

This appendix is focus on a derivation of the first expression of partition function ^ based Wang's formal framework 

Q. 

Note that the Hi matrix is able to be taken as off-diagonal based on Wang's proof 1], from Eq.Q it follows 

71,72 1=1 ^ 

It implies that the contraction and anti-contraction of off-diagonal elements skill using in [l[ is largely simplified when 
Hi matrix is taken as an off-diagonal form. 

Since the existing the factor S-yj^S^j^^-^ in the A^'^ expression ([3]), it leads to that E^^ = E^^^-^ = E^ after summation. 
It implies that there is, at least, an obvious singular point among the beginning and ending terms in the summation 
over i from 1 and / + 1 for a fixed I > 1. To remove it, we can combine this two terms, and introduce an infinite small 
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ie to represent the difference between E-^-^ and -^7,^1, and then replace E^-^ by E^, again finally set e ^ 0. That is: 



/+1 ; 

k J=2 j = l 



lim e"'^-'' { 
£^0 



dz 



.i=2 



(i£) (-i£) 



771^771 + 1 



(-B^ - £^7, - z) 



(A2) 



2=0 



The last step has directly used the definition of derivative. 

Furthermore, let us first consider how to deal with the terms in the summation over i from i = 2 to i = Z of Eq.Q. 
By summing the index 71 and ^i+i we have 



7,7i,-,7, + l-2 -Q {^E^^_E. 



-5 5 TT o^^'^'+i = 'S^ 
"771 "71+17 1 1 y / ^ 



(£'72 - 



^772^727^;, 



7j; 



- E^Y{E^^ - E^^) (£'73 - £7)2 (£^3 - £;^J 



^772^7273^737^^3 



7.72,--- ,7; «=2 













nV7. 








_J=2 







^772 



7-1 
.1=2 



7i7i + l 



5^'^ (A3) 



where 9{x) = 1 if a; > 0, otherwise 9{x) = 0. 

The simplest case is that 1 = 2. Only there is the first term in above equation (|A3|) . Interchanging the dummy 
index 7 ^ 72, we have 



E 



(^^72 - E^r 



gill gill _ ' Q-i-E-,t 



E' 



\ 

dz \ {E12 ~ E^ ~ z) 



9' 



7,72 ^ ' ' 7,72 

Similarly, for ? = 3, we set {72,73,7} {7,72,73} for i = 2, and {73,7,72} ^ {7,72,73} for i = 3. Thus 



(A4) 



-iE^.t 



E E — 

7,71,-, 74 .=2 -Q 



"'^771^74 7 Jl 5 



7i7i+l 



1]J 



E 

7,72,73 

= E 

7,72,' 

= E 



1 1 



gill glll-i glZl 



[Ey - E^2Y{E^ - £73 ) {E^ - Ery^ ) (E^ - E^^Y 



gll2 gl2l3 gl3l 



^-lE^t I gi2t 



7,72,73 



_d 

d^ 



{E^ - Ej2 - z)iEi - £73 - z) 



gll2 gl2l3 gl3l 



Z=0 



(A5) 
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The skill of dummy index transformations can be continuously used when I > 4. For i = 2 and to i = 
they are respectively {72, 73, ' ' ' , 7i, 7} ^ {7, 72, • • • , 7i-i, 7h } and {7;, 7, 72, • • • , 7i-i} {7, 72, • • • , 7;}- For 
the other i {I — I > i > 3), our dummy index transformations are taken as {7^,7^4-1,- •• ,7/,7,72,'' - ,7i-i} 
{7,72,' •• , 7;-i+i, 7;-i+2, 7/-i+3, • • • ,7;}- It is easy to prove that, under above index transformations, 



gii'2gi2i3 . . . gin invarlaut in form. Obviously, i = 2 and i = I terms are transformed as following form 



k=3 



Lfc=2 



iE^-E^,r 



[E^, -E^YW [E^, - E^, ) {E^ - E,,Y W [E^ - E^ 



k=2 k=3 

While the other i from 3 to / — 1 > 4), our dummy index transformations lead to 



i-i 

E 

1=3 



n (^7.-^7.) 



l-l 

E 

i=3 



{E,~E,,f 



n (^7-^7,) 



Thus, when I > 4 



(A6) 



(A7) 



(A8) 



E E E- 

i=4 7.7l-7,H-l^=2 -Q [E^^-E,;^ 



.fe=i 



ikik+i 



'^771 '^771 + 1 



l=A 7,7273 ■■ -7! 



X{{E^-E 

J=2 



{E,-E,,r 



+ E' 

i=3 



+ - 



I 



{E^~E,,yl[{E^~E,J 

d 



E E 

l=i 7,7273---7( 



k=3 



l-l 

1=2 



e 



dz 



n 

i=2 



,7i7i+l 



1 



(-E'7 ^ ^7j ^ ^) 



5 



772 



z=0 



7i7i+i 



lis"" 

.1=2 



^7^7 



Substituting eas. (jAllA2l A4IA5I A9p into the expression of partition function, we obtain 

772 7273 . . . „7!7 



^^77 = 5^0--^* + Y: e --9' 



d 

d^ 



1=2 1,12, ■■■ ,7i 
/ 



3 ' * J] 



1 



J=2 



- i?^^. - z) 



n 



1 



J=2 



{E, - E,^ - z) 



E' 



-\E,,t 



1+E E (-w^^s'^^^^-'-s'^'^ 



;=i 7i,--- ,71 



1 

n(^^7-£^7.) 



{E^-E,J^ 


I 

n (^^7-^7,) 






'l-l 





z=0 



(A9) 



(AlO) 



in which, we have reset the summation indexes. 



10 



APPENDIX B 



In this section, we would like to further remove all apparent singular points in the expression of partition function. 

It is clear that there are still the singular points in the expression (|A10|) although they are fake. In fact, these 
singular points as well as the relevant terms in form are similar to those have been removed in Appendix |^ except 
for the delta function factors. Consequently, to further deal with these unexpected singular points, we need to pick 
out the terms when = E^. . Its method is just like the things that the ref. [i| has ever done, that is, we start from 
the identity in the sense of summation: 



1 — ^^'yi ^TTj 



(Bl) 



and rewrite a summation as 



^ ^ / ['''71 '''73] ^ ^ / [-^7 7 ■^-fi] ^■yfi ~l~ ^ ^ / [■^-fj •^7s] ^77i 



(B2) 



The first summation of right side of above equation is picked out, in which, x-y = x-y-. Extending this method to a 
Z-fold summation, we chose m elements from a set £ = {1, 2, • • • , 1} to form a subset {p'}m, so that 



1 - H 51 V}-7 
m=0 {pi}„ 



'{p'}n.7 



(B3) 



(B4) 



where {q^},n — ~ {p'}m- Obviously 



, pj„ e £ and p- < if i < j} 



G (£ - {p'}m) and q\ < if « < j} 



(B5) 
(B6) 



In fact, if p\ — = 1, then g"'^'^^ 5^.,y5^^. = since g'^^ = has been taken here. This result leads to that the number of 
subset {p'}m element with contribution is not larger than [{I + l)/2]. However, in form, we can keep these vanishing 
terms. 

/ 

In order to pick out the singular terms from the summation e~'^"'*]^(i?^ — E~f.)~^ we product it by '^{p!},^^ 

7i>---:7l i=l 

so that there are, at least, m obvious singular points within it since <5{p'}„7 contains m delta functions. Obviously, 
TO = case is not needed to considered since there is no the singular point. 

Similar to the skill used in Appendix |^ we rewrite it as (to + 1) terms, in which, the first term does not involves 
the delta function action, but each in the other to terms respectively absorbs the contribution of each delta function, 
that is 



1 



{P'}m7 



TO + 1 



I I 



'{p'}™7 



(B7) 



In order to remove the fake singularity, we introduce m infinite small numbers £m (m = 1, 2, ■ • ■ , so that E^ — E^^^ 
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iSm- Using the fact that E^^, — E^^, = —i{si — ej) and noting the action within S^piy^^ , we have 



HiE.-E,^) 

.3=1 



Q — iEyt 

hrtn -v = T lim 

£2-^0 
Sm^O 



l — m 



ir-Sk n - ^k) [] {E^ - E^^, - i£fe) 



'{P'}m7 



fe=0 



m + : 



(B8) 



(B9) 



where we have denoted the gj- G {q''}m [j = 1)2, •• • ,Z — m), rewrite the summation or production over the index 
belonging to {q^}m (or ^ {p'}m) as from to ^j.^, and also define B}. by 



Bk 



rife) \{{E^-E ) 
i=i j=i 

1 



l — m 



(fc>l) 



(BIO) 



(Bll) 



n (^i - sfe) n " -^-r,' " ^^fc) 



Now let us prove the following equation 



11(^7-^^7.) 

.3=1 



'{P'}m7 



(-1)™ d" 



(m + 1)! d^" 



(^^^^^7^ 



'{P'}m7- 



(B12) 



z=0 



Firstly, we calculate the limitation ei — > 0. Only ei is an obvious singular point in Bq and Bi, that is 



}imiJ2B^~S^,>y^, = Hm 



-eit 



ei n n - E^,^. ) ^1 n - n - ^^7,, - 



i=2 



3 = 1 



i=2 



l — m 



y 

The last summation term appears only when m > 2. Actually, by setting 



3 = 1 



'{P'}m7 



'{p'}m7 



(B13) 



/L(^) = 



1 



fc — III 

Wie^-x) Yl{E,-E,^, -ix) 



(B14) 



next doing its Taylor expansion 

ftM = flm{0) + 
1 



x=0 



Ci2 ^ 



x=0 



i=k j=i ' 



1^ 

2! dx2 



d 
d.T 



/ — m 



\{{e:-x) \[{E-,-E,^, -ix) 

i—k j—i 



and then substituting k — 2 result into Ea. (|B13| ). we arrive at 



lim f V B, 



-El* 



i=2 



£1=0 



fc=2 



i—2.i^k J — 1 ^ 



'{p'}™7 



From the definition of derivative of a function, this result is very obvious. 
Secondly, let us find the limitation 62 ^ 0. It is easy to see 



{P'},T17 



\k=0 



lim 

£2^0 



'{P'}m7 



fc=3 



Hp'}. 



where fi„^{x) is defined by Eq. (|B14p . From 



(-i)d^ Li=o = -^/L(o) - ^ 



dei 



£1=0 
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and Eq. (|B15p but taking fc = 3, it follows that 



\fc=0 / 



(-1)2 



£2* 



i=3 



'{p'}™7 



fc=3 



l—m 



£2=0 



<5{p.}„^ (B19) 



In particular, 



\fe=0 / 



4)2 d2 



2! dei 



-£2* 



1-2 



{p'}27 



(B20) 



£2=0 



Analyzing the above derivation, we have seen the fact that the singular points are able be removed by finding 
limitations step by step. Every limitation calculation increases one order of derivative in the first term and decreases 
a term in the last summation. When all limitations are calculated, this expression becomes a pure m order derivative. 
Without loss of generality, for n < {m ~ 2) we assume that 



£2^0 \fc=o y 



{P'}»7l7 



(-1)" d' 
n! de 



,=0^{p'}m7 



£„-»0 



+ 



E 

k=n+l 



-£fct 



l — m 



i—n-\-l,i^k J — 1 ^ 



(B21) 



we can have 

i™0 (E^MV}™7 



£2^0 \fc = 
£„+1^0 



lim , 

£„+i-*o n! de: 



E 

k=n+2 (•_g^_'Jn+2 



£„=0 



"■n+1 J 



-£>:* 



2— n+2,2^fc 



Using the fact that 
d^ 



de: 



„ [fhn i'^n)] Ig^^Q 



d" 

de" 



1 



1 dJ 
j! dT^ 



flm (^") 



i=i 



^ n! 1 _d^ 
Z^TTTh^ dZ 



<5|p,}„^ (B22) 



£n=0 j=o ^n+1 



flm (^") 



£,1=0 



E ,1 A J ■^i'^^^'^") 
J=0 



-ri+1 



(B23) 
(B24) 



£„=0 
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we have 



lim V 



Bk (^{p'},, 



-1) 



ri+l 



(n + 1)! 



jn+l 

+ lJlm \^n+l> 



de 



n+1 



'{p'}™7 



£,1+1=0 



1 6"^'='' 



J{p,}„^ (B25) 



Finally, set n = m — 1 and then use the same method we can finish the proof of Eq. (|B12[) . Obviously 



£2^0 \fe = / 
£m->0 



= lim 



(-1)^ 



(m - 1)! 



jrn— 1 



de^ 



-1 flmi^rn-l) 



-m l-r. 



£to_1=0 



'{P'}n.7 



(-1) 



d'^ 



l — m 

1[{E^~E -is) 



(B26) 



£=0 



Substituting it into Eq. (jB9p and setting ie — z follow the conclusion Eq. ()B12[ ) 

When the set {q^}i is not an empty set, the product or summation over this set is well-defined. Usually, in the 
formal expressions, the product over an empty is thought of 1 and the summation over an empty is thought of 0. At 
the above sense, m still can take But g~''^ = has been taken, thus 



9"'g' 



■9^'''hp'}n =^ 



(B27) 



Finally, we arrive at 



5: A- . ^e- 



oo l — l 

E E E E (-^)°" 

1 = 1 m=0 77172---7I 



(-1)" d™ 



n (E, - E^^ - z) 



z=0 



g77lg7l7....^7,7^^^,^^^ 

It must be emphasized that all obvious singular points have been removed in above form. 



(B28) 



APPENDIX C 



Now, our task is to continue to derive out the final expression of partition function that is written as a series of 
power of a kernal function as well as its derivative. 

It is clear that is a subset of £ = {1, 2, • ■ • , Z}, and has, at most, / — 1 elements. Its every element can be 

taken as an end point so that C is divided into m + 1 subsets = {Pn~i + 1>pL-i + 2, • • • .p^} [n = 1, 2, • ■ • , m + 1) 
except for pm = I case. But this exceptional case does not really appear since g'^'^"^^5^pi]^^ = 5''^'^{p'}„7 = by 
using the fact that g'^^ has been taken as off diagonal, but in form, we simply set C^+i to be an empty set. 

Obviously, the cardinal numbers of are respectively In = Pn — p\i-i: in which Pq = and pj„+i = I are defined. 
For unify the denotation we need to define the subsets that do not contain the p\ , they are — Ln ~ Pn and in the 
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exceptional case when pm = I, C-'m+i = ^- Their cardinal numbers are respectively Z'^ = ^„ — 1. But in the exceptional 
case when p„ = I, we still set I'^j^i = 0. 

After summing all delta functions within <5-y^^,j ^ and grouping the relevant terms, we have 



m+1 



n 



[E^ - E^^ - z) 



E 



7„!. 



n 



(i;^ - - z) 



(Cl) 
(C2) 



It must be emphasized that = will lead to appear the factor g ^^+^5^ ,jSj , j = g'''^Sj ,jSj , 7 = 0, and no 

any — ^ is grouped into the square bracket. Note that the dummy indexes can be changed according to 

the given rules, we define 



7l,--- ,7i' 



L {E^ - E^^ - z) 



(C3) 



and then 



I 



77!.+ 1 



771 n7l72 . . . „7i7 



n 



7l.--->7i i=l,i^{p'}^ 

Again substituting it into Eq. (|B28p . we arrive at 

oo I— I 

^^77 = eh^)^ 

7 7 i=lm=0{p'},^ 



(C4) 



(-1) 



(777+1)! dz"- 



n=l 



'm+1 



(C5) 



2=0 



The summations over all set {p'}m can be changed into the summations over l'i,l'2, ■ ■ ■ , l'^, but there is a limitation 
that X]^=i ^" — ^- Interchanging the summations over I and 777, then / begins at 777 + 1, setting I'^j^i = I — m, and 
noting I'jy^^i can be taken up to infinity, every Z'^ (77 = 1, 2, • • • , 777) also can be taken up to infinity. Therefore 



Y^A^^ ^Y. + E E E (-iOe-'^^* 

m— ^ ^ 

where we define the function Rj{z) by 



(-1) 



(777 + 1)! dz 
i?7+i(z)e'^*] 



m+1 



.n=l 



z=0 



z=0 , 



(C6) 



i?7(.) = E4''w = E E 



g11lglll2 . . . gin 



1=1 



1=1 7172 -7, #7 '■^'^ ~ ^"^1 ~ ^'''-^'^ - -6^72 - ^) • • • (-5^7 - ^71 " ^) 



It is a kernal function in our approach and expressions. 



APPENDIX D 



Furthermore, we can expand the partition function as the power series of time t. 



(C7) 



E-yt 



1 + E 



i-ity 



n=l 



(Dl) 
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where the coefRcient of {—it)^ is as below: 



C2 = E 



ru— 71+1 



^ (m + l)(m - n + 1)! 

— n— 1 ^ ^ ^ ^ 



jm— n+1 



(D2) 



z=0 



It is easy to prove by expanding e'^* and noting the derivative at z = 0. 
From the definition of partition function, we also have 



7 7 



7 7 

Here, we have rewritten the Hamiltonian eigenvalues as 



After finishing the proof of the following relation 



we can obtain 



Obviously 



(D3) 

(D4) 
(D5) 
(D6) 



^7_ V (-ir 



^ (m+l)! 

m— 



dz™ 



(D7) 



2=0 



Now let us begin our task. It is clear that if C^^_^i = (C'Z) ^ ^'^ verified then Eq. (|D5[) is proved. From the 
definition of CJJ, we see that. 



CZxCl = \ E 



^-^ (mi + l)(mi — n + 1)! 



jmi— n+1 



E 



, ^„ (to2 + 1)! 

M-71 

E (-1)"-" E 



dz™2 



^-'^ (M - m){M -m- n)\{m + 1)! 



" jAf — n — m 



^2 Af — n — m 



dz™ 



(D8) 



2 = 



The last step has used the index transformations mi + m2 + 1 and m2 — > m. While, we also know that 



ci^i = E (-1) 



Af-r 



(n+1) 



(M + l)(M-n)! 



]A/-r; 



dz 



Af- 



(D9) 



2=0 



So, to prove Eq. (jD5|l . we should prove the equation below (where n > 1). 



jAf-n 



M-n 



n M + 1 



(M-n)! 



dzA'-f 

Or for simplicity, set M — n = k 



^ n + 1 M - TO (m + l)!(Af - n - m)! \dz*^-"-™ 

m=0 ^ 



3M—n—m 



-R 



M-r 



R':;+\z) (Dio) 



d^ 



-R. 



fc+n+l 



w-E 



n fc+n+l 



fc! 



m=0 



71 + 1 fc + n — TO (to + l)!(fc — m)! ydz'' 



(^) 



dz"^ 



dz™ 



(Dll) 
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Note that it is not Leibnitz' rule for the nth derivative of a product since the differential functions connect with the 
summation index. It is clear that when k = 1, we can directly verify Ea. (|Dlip is valid. That is 



E 



n 1 + n + 1 



^ n + 1 1 + n - TO (m + 1)! (1 - to)! V dz^" 

m— ^ 

= {n + 2)i?"+i(z)^i?^(z) = 



-R 



1+n— rn 



D 

dz" 



m+l 



(^) 



(D12) 



Assume up to a given I Eq. (|Dlip is correct, and then, in the case of Z + 1, we can also prove it is correct as below. 
After the one order differential, using the Leibnitz' rule for high order derivative of a product, it follows that 

I 



d'+i 
dz'+i 



n 



i+i 



dz 



(D13) 



Since Eq. (|Dll[ ) is assumed to be correct up to a given I, Eq. (jD10[) is valid up to Af = ^ + n (rt > 1). Thus we can 
make use of Eq. ljDlOp to replace the {I — j)th derivative in above equation except for j = 0. 



dz'+i^ 



; + l+ri+l 



1-3 

E 



i + i 



(i-jV- 



-R 



( d" 



-R 



Vi + 1/ \l-mj (to + - j - my. 
d^ 



m+l 



Interchanging the summations over j and to, we have 



dzJ V^"dz^^ 



d' 



^"dz-^^ 



^-'o ('-"i) (™ + l)! 



E 



d^ 
dzJ 



dz'^ 



-R. 



l + n + 2 / d 



n+1 



dz 



Specially, the terms in the square bracket can be simplified as 



• (D14) 



(D15) 



E 



1 



R 

jy Vdz'--?-™ ^ 



l — m 



^ ( i?"+A 



dz3 



dz 



E 



^ j!(;-TO-j)! Vdz'--'-™ 



-Rt 



l — m / ^ 

^ VTTT 



1 



dJ 

dzj" 



R"—R^ 

' dz 



TDTl p 

j!(;-TO-j)! Vdz'-J-™^''' ; dzJ V ''dz 



-i?: 



d^ 



(/-to)! 



d'+i 



(/ + 1 + - to) dz'+i- 



R 



/+l+n— m^ 



1 / d' 



l-r. 



+ 1 ^ (7 + 



+ 1 ^ (i + 1) iK^ -m~ jy Vdz'--?-™ 



-r: 



I — m 



n+1 \dz^ 



— R^ 



n+l 



1 



(/-to)! 
1 



d'+ 



l — m 



(/ + 1 + n - to) dz'+i-™ ^^'^ 



/+1+71 — ^7^^ 



dzJ+i 
1 / d'-™ 



n + 1 Vdz'-™ 



1 — 771 



1 



(n + 1) (/ + 1 - to)! 



-R. 



/+l+n— m 



(/ + 1 - to) 

1 



dz'-' 



-R 



d 



/ j; + l-m 



dz^ 



n 



n + 1 (/ - to)! (/ + 1 - to)(/ + 1 + n - to) dz'+i-™ ^ 



1 



1 



/ d'+i- 



R 

+ 1 (/ + 1 - to)! Vdz'+i-™ 



_^ pn+1 

dz^ 



p 



,n+l 



(D16) 
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Substituting it into Eq. (|D15p . we immediately arrive at 



d'+\^,+i+„+i ^ ^ n ; + l + n+l 



(/ + !)! 



dz'+i ^ 



'^n + I I + I + n - m [m + + 1 - m)! Vdz'+i" 



-R 



/ + l+n — rn 



(^) 



E 



1 l+n+2 



(/ + !)! 



n+1 l-m (to + 1)!(? + 1 - m)! Vdz'+i-™ ^ 
d 



m=0 

; + n + 2 / d' 



'7 



-i?: 



n+l 



n + 1 \dz' y dz 

E'^""" n / + 1 + n + 1 
n+U + l + n- TO(TO + + 1 - to)! V dz'+i-™ """^ 

m— ^ 



(/ + !)! 



/ d'+i- 



; + ?i + 2 
n + l 
i-i 

-E 



1 



1 + 2 Vdz'+i y ^ n + l Vdz'"'' 
(/ + 1)! /d'+i 



/ d" 



m=0 



-i? 

l-m{m+ + 1 - to)! Vdz'+i"™ 7 



d"^ 



1 — ^T^' 

dz™ ^ 



Now, the problem is changed into the proof that the below expression is zero. 



1 



d'+i 



1 + 2 Vdz'+i ^ 
/-I 



7 

(/ + 1)! 



n-l)[ --fRl+ I ^i?: 



?ri+l 



m=0 



R 

Z - TO (to + 1)!(? + 1 - to)! Vdz'+i-" 7 



dz' ^ y dz ^ 
d'+i 



dz' 



-R. 



m+l 



Vdz" ^ 



(D17) 



(D18) 



Similar to deal with Eq. ljDlOp . that is, carrying out derivation like Ea. (|D17[) . we can obtain the {I + l)th the derivative 
in the first term. 



(/ + 1)! 



d'+i- 



^ Z - TO (to + 1)!(? + 1 - to)! Vdz'+i" 

771 — ^ ^ ^ ^ ^ 



— R^ 

-m 7 



\dz 



(D19) 



Substitute it into Ea. (|D18| ). we immediately see the expression (|D18[) is zero. It means that (jPlip is proved to be 
correct for Z + 1. Therefore, we have proved Eq. (jD5[) 
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